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ABSTRACT 


We have detected and spectrally resolved the mid-IR hydrogen recombination lines H6a 
(12.372pm), H7a (19.062pm), H7p(l 1.309pm) and H8y(12.385) pm from the star MWC349A 
This object has strong hydrogen maser emission (reported in the millimeter and submillimeter 
hydrogen recombination lines from H36a to H21a) and laser emission (reported in the H15a, 
H12a and H10a lines). The lasers/masers are thought to arise predominantly in a Keplerian disk 
around the star. The mid-IR lines do not show evident signs of lasing, and can be well modeled 
as arising from the strong stellar wind, with a component arising from a quasi-static atmosphere 
around the disk, similar to what is hypothesized for the near IR (<4pm ) recombination lines. 
Since populations inversions in the levels producing these mid-IR transitions are expected at 
densities up to ~ 10" cm' 3 , these results imply either that the disk does not contain high-density 
ionized gas over long enough path lengths to produce a gain ~1, and/or that any laser emission 
from such regions is small compared to the spontaneous background emission from the rest of 
the source as observed with a large beam. The results reinforce the interpretation of the far-IR 
lines as true lasers. 


Subject Headings: circumstellar matter - stars: emission-line. Be - stars: individual (MWC349A) - 
stars: mass loss - infrared: sources - infrared: spectra, masers 
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I. INTRODUCTION 


MWC349A, the brightest radio continuum star in the sky, has a strong bipolar wind with 
mass loss estimated at ~10' 5 M 0 yr ’ (Cohen et al., 1985), or even as much as ~ 9xl0' 5 M 0 yr'' 
(Martin-Pintado, Bachiller, Thum and Walmsley, 1989). It is surrounded by a Keplerian disk of 
ionized gas seen edge-on (e.g., Hamann and Simon, 1986; Gordon, 1992; Thum, Martin-Pintado 
and Bachiller, 1992; Ponomarev, Smith and Strelnitski, 1994, hereafter Paper I; Strelnitski, Smith 
and Ponomarev, 1996, hereafter Paper III). The source is so far the only confirmed emitter of 
strong hydrogen recombination line masers and lasers in the millimeter (e.g., Martin-Pintado et 
al., 1989a, b) and far infrared (Strelnitski et al., 1996). Laser emission (for simplicity we will use 
the term laser when both maser and laser are meant) is seen in a series of Hna lines (that is, 
transitions with n+1- n) starting from about the H40a line at 135GHz (2.22 mm) and extending 
down at least as far as the H15oc line at 169pm (Strelnitski et al., 1996, hereafter Paper IV), with 
likely amplification seen even at the HlOa line at 52.5pm (op cit.), and in one Hnpline, the H32P 
at 366 GHz (Thum et al., 1995). 

Continuing from both ends of this string of lasing lines are many other recombination lines 
which have been reported from this source. At radio wavelengths lines have been reported out to 
the H92a line at 3 cm (Rodriguez and Bastian, 1994), and from the opposite extreme they 
continue through the near infrared H4a (4.05pm, Hamann and Simon 1986) into the visible 
H2a=H a (6524A, Kelly, Rieke and Campbell, 1994). In addition in the near IR higher order 
recombination lines are seen -- Hnp, Hny, and above (e.g., Kelly, Rieke and Campbell; Hamann 
and Simon, 1988). Absent from this long series of transitions, covering over four orders of 
magnitude in wavelength, are the lines in the critical region between 4pm and 52.5pm. The mid- 
IR range is important for two reasons. First, somewhere in this region the detected lasing 
disappears, despite theoretical indications that the populations of the corresponding levels are 
inverted. Observationally the strength of the lasing is seen to steadily decrease to the short 
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wavelength side of the peak in the submillimeter where the excess above the expected 
spontaneous emission is a factor of ~200 (Paper III and Paper IV) to small values (~x6) at 52pm, 
just as on the long wavelength side the maser action appears to die off longward of about the 
H40a line (e.g., Paper III). The second reason for its importance is that the mid-IR lines, if 
lasing, would sample higher density gas thought to arise at the inner edges of the disk, closer in 
than the ~5AU modeled for the HI 5a line ( Paper IV). In this paper we report on velocity 
resolved infrared spectra of four recombination lines in this interval, and discuss their relative 
intensities and profiles. 

A great advantage of the masing transitions has been that the strong, velocity resolved 
emission of this sequence (albeit complicated by the fact they show substantial time variations) 
enables the detailed modeling of the environment. The bipolar wind is recognized as being the 
source of the strong radio continuum flux, and of the non-masing centimeter radio recombination 
lines (Rodriguez and Bastian, 1995). Although early suggestions also proposed it as the source of 
the intense maser lines (Martin-Pintado et al., 1989a,b) it is now clear that at least the strong, 
double-peaked emission arises from the disk and not the wind (e.g., Gordon, 1992; Thum, 
Martin-Pintado and Bachiller, 1992). A portion of the "pedestal component" may also arise in the 
wind. A specific model has been proposed for such "photoevaporated" Keplerian disk structures 
by Hollenbach et al. (1994), and has been used to obtain physical parameters of the disk (Paper 
IV). The non-masing near-IR and optical hydrogen recombination lines, as well as a plethora of 
other atomic lines, have been attributed to various zones of ionization along the radius of the 
Keplerian disk (Hamann and Simon, 1986), as have the lasing far-IR lines (Paper IV). 

Maser and laser amplification in the submillimeter and FIR lines implies quite high 
densities, N e > 10 6 ' 7 cm' 3 (e.g.. Paper I; Strelnitski, Ponomarev and Smith, 1966, hereafter Paper 
II). The presence of very high densities of circumstellar ionized gas in this source, even as high as 
N e .slO 10 cm' 3 , was proposed by Hamann and Simon (1986; 1988) to explain the emission from 
[Fell], [Felll]. Kelly, Rieke and Campbell (1994) also found evidence for densities N e slO 6 cm' 3 
in their near-IR spectra. Earlier, Smith, Larson and Fink (1979) suggested lasing in the 
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abnormally bright infrared hydrogen Pfp line from Orion/BN, but noted that at the high densities 
implied (>10 10 cm' 3 ) improved models for the hydrogen levels populations were needed to model 
the line strengths. Hummer and Storey (1992) and Storey and Hummer (1995) have recently 
provided such calculations, with densities that extend beyond N e a 10 10 cm' 3 . We have used their 
tables (Papers I, III and III) to model the maser and laser emission from MWC349, and in 
particular to confirm the prediction that population inversions persist to these high densities. As 
the transitions having maximal laser gain shift steadily towards lower quantum number n as the 
density increases, at densities above Ne > 10 7 ' 8 cm' 3 the far and mid infrared Hna lines are 
expected to be efficient amplifiers. Moreover the very large gain coefficients expected from even 
modest sized regions with such high volume densities. Since the gain is nominally proportional to 
the density squared, the H7a line for example has its characteristic length as only 3xl0 9 cm at 
N e alO 10 cm' 3 (c.f. Paper II), suggesting these lines might be detectable and useful probes of the 
inner few AU. 

We observed a set of mid-IR recombination lines to determine whether, given the 
theoretical possibility of inverted populations and high gains, the trend of observable lasing lines 
continues down into the mid IR in MWC349A, and to determine whether this enables us to refine 
further the parameters of the Keplerian disk. The result might help shed further light on the 
remarkable fact that despite the ubiquitous inversions predicted in hydrogen only MWC349A to 
date shows unambiguous and strong laser emission. In this paper we report our results: velocity 
resolved observations of four mid-IR hydrogen recombination lines from MWC 349A — the H7a 
line at 19.062pm, the H6aline at 12.372pm, the H7pline at 11.309pm, and the H8y line at 
12.385pm. 
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II. OBSERVATIONS AND ANALYSIS 
a) Instrumental and Observational Summary 


The hydrogen recombination lines were measured from the NASA Infrared Telescope 
Facility on the nights of June 7 and 8, 1994, with a followup search for time variability one year 
later, in July 10, 1995. Table 1 shows a log of the observations. The Irshell spectrometer used 
either a 1" or 2" slit width (depending on seeing conditions) and had a resolution of 20-30 kins' 1 . 
Flatfielding was done using a dome-temperature card chopped against the sky. Atmospheric 
corrections were done by observing the standard x Cyg immediately before and after each 
observation, flatfielding in the same way as the source spectrum, normalizing the average 
spectrum to unity, and dividing it into the source spectrum. Flux calibration was done by 
normalizing the continuum flux of the observed spectra to an average spectrum generated by 
fitting the photometric observations between 9 and 20 pm of Simon (1974), Simon and Dyck 
(1977), Herzog, Gehrz and Hackwell (1980), and Aitken etal, (1990). Velocity calibration was 
done using the telluric lines in each spectrum. In general the fluxes are precise to 5-10%. Data 
files were processed using Irshell's own software; further information on the instrument can be 
found in Lacy et al. (1989). 

Figure la shows the 19pm spectrum, as well as four gaussian profiles into which it can be 
decomposed. Figure lb shows the 12.4pm spectrum containing both the H6a and H8y lines, 
along with its gaussian decomposition. Figure lc shows the 1 1.3pm spectrum containing the 
H7P line. The locations of weak telluric features are also marked throughout. Table 1 lists the 
observed parameters of the four measured lines and their decompositions; repeat observations in 
1995 showed no variations from these. 
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b) Lasing, or Optically Thin Emission ? 


In presenting their analysis of the FIR lasing lines H15a, H12a and HlOa, Strelnitski et 
al. (Paper IV) point out the trend, first described by them in Paper III, that the integrated Hna 
lines fluxes in the millimeter and far IR regime (H36a - HlOa) deviate in intensity from that 
expected from optically thin, spontaneous emission. The latter is shown to be a smooth curve 
nominally proportional to n' 6 to n' 8 , where n is the principal quantum number, and the range in 
exponent is due to the increasing importance of free-free absorption at the longer wavelengths. 
Strelnitski et al. (1996) use this excess emission over the optically thin value to show that the FIR 
lines measured from the KAO are lasing, but with a smaller amplification than the maximal 
amplification values seen in the submillimeter. The integrated values of the two mid-IR a -lines 
presented here fall close to the expected spontaneous emission values predicted by the 
aforementioned dependence, providing evidence that the net emission in these lines has indeed not 
been significantly amplified. In the next section we show this in more detail. 

c) Emission from the Wind and the Disk 

The source structure of MWC349A is complex. Hamann and Simon (1988) delineate at 
least five different emission regions, and Kelly, Rieke and Campbell (1994) also find evidence for 
each of them in the near-IR. Rodriguez and Bastian (1994) observed the H92a line and conclude 
the bipolar wind is itself rotating. Finally the variability and asymmetries of the maser line profiles 
(e.g., Thum et al., 1994) indicate both complex substructure in the disk and possible infall. Many 
of these structural subtleties are still under debate, but we think it is safe to admit that at least two 
distinct regions contribute to the recombination line fluxes we present here: the bipolar wind, and 
the quasi-static, ionized atmosphere of the Keplerian disk. 

In order to estimate the contribution to the mid-IR fluxes from the bipolar wind we 
normalize the calculation to the flux in the H92a line measured by Rodriguez and Bastian (1994), 
and Rodriguez (private communication) at (7.6 ± .8) xlO' 27 W/cm 2 . Unlike the H66a line, for 


7 



which there is at least one suggestion it has some component of masing in the wind ( Martin- 
Pintado et al., 1993), the H92a line shows no evidence for anything except normal emission from 
a strong, ionized wind. Indeed the radio continuum emission from MWC349A follows an L v « 
v 0 6 law down to wavelength of about 100pm, as is characteristic of such an optically thick flow. 
The hydrogen lines from such a wind sample surfaces of optical depth x« 1. Their (normalized) 
fluxes are calculated assuming the wind has a constant velocity flow V=V i (r/r i ) s with s=0 (e.g., 
Smith, Fischer, Geballe and Schwartz, 1987; Panagia and Felli, 1975). The results are plotted as 
the solid line in Figure 2; fluxes for the entire series of Ha lines from H4a to H 16a (wavelengths 
indicated by triangles on the x-axis), including those lines not yet reported, have also been 
calculated and plotted. The measured fluxes are remarkably close to those extrapolated from the 
H92a flux, but they are stronger by a factor of 1 .8 for H6a and 2.0 for the H7 a lines. The same 
model finds the H7p and H8y line fluxes are consistent with a wind to within their measurement 
errors, and these results are shown in Figure 2 as small open triangles that virtually overlap the 
measured data points (the curve refers only to the a-line transitions). The Bra line flux is also 
about 90% higher than the wind prediction. The lasing H15a line is, of course, 35 times brighter 
than this model. The 10% uncertainty in the H92a flux we have taken as a reference is not 
enough to explain these deviations from the wind prediction. 

There are several possible sources for these deviations from wind model predictions. 
First, even a small change in the wind parameters, for example a modest acceleration like s=0.5 
instead of s=0, would be more than enough to boost the near and mid-IR line fluxes when 
extrapolated from the H92a line which is produced so far away (see the discussion and relations 
in Smith, Fischer, Geballe and Schwartz, 1987). However, Hamann and Simon (1988) note that 
their [Fell] emission might suggest deceleration : this would push the predicted fluxes of the H 
lines even further from those observed. Alternatively, as noted above, the wind itself may have 
very complex structures and/or motions. We have generally adopted (see also Strelnitski et al., 
1996) the generic "weak wind" model of Hollenbach et al., (1994) for this source to describe 
aspects of the disk flaring and wind, but note that they too call attention to the fact that for 
MWC349A their model has difficulty describing the observed emission within the inner 30 AU. 
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We note that their model predicts an n<*r' 3/2 dependence to the flow, whereas the continuum data 
seem to be consistent with a l/i^wind (e.g., Cohen et al.) which we have adopted to fit the line 
fluxes. 


We think that the most probable cause of the deviations from wind models is that the 
morphology of the emitting ionized gas around MWC 349A is not limited to the wind. In 
particular, in the Hollenbach et al. model the stellar ionizing photons maintain an HH-region like 
atmosphere above the disk, including the region of the inner zone where the near and mid-IR lines 
are expected to arise. Hamann and Simon (1986) attribute the location of the Bra line to such a 
consistent location: the double-peaked profiles and/or asymmetries of the shorter wavelength lines 
are strong indications of their originating in a rotating disk rather than in a bipolar wind. We 
therefore modeled the residual (i.e., non-wind produced) fluxes as arising from a static, optically 
thin HII region, and calculated the values from the Hummer and Storey Case B emissivities for 
the case of T=10 4 K and N e =lxl0 7 cm' 3 ; the conclusions are not particularly sensitive to the exact 
values of either T or N e . Figure 2 plots the sum of the contributions to the lines, wind plus "HII 
region", where we have normalized the latter contributions by the total flux in the Bra line minus 
the predicted wind contribution. 

Using the canonical value for extinction in this source, A v =10 mag (Cohen et al), the 
total Bra flux is 3.0xl0 17 W/cm 2 (Hamann and Simon, 1986), with about half the flux contributed 
from the strong wind. The predicted "HII region" contributions to the mid-IR lines, when 
normalized to the remaining Bra flux, are rather close to those observed: only about 20% above 
the observed values for both the H6a and H7a lines. If the assumed extinction is too high, the 
intrinsic flux of Bra is reduced and all the predictions fall even closer to the measure values. 

The current extinction value is derived by Cohen et al ' s analysis of the photospheres of both 
MWC349A and MWC349B. They find the extinction to MWC349A is about 2 mag larger than 
that for MWC 349B, and they speculate that this extra extinction is likely due to the effects of the 
edge-on disk, while that of the other 8 mag of extinction is probably foreground dust. If the Bra 
and mid-IR lines arise in a region outside of the disk, they might not experience all of this 
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additional extinction. Given the other general uncertainties in the physical structure of this inner 
region, and the wavelength dependence of the extinction itself, we conclude that an extinction 
between A v =6-10 mag can adequately describe all the mid-IR line fluxes. We note that Kelly et 
al.{ 1994) also estimate an A v ~ 1 1 mag from the near-IR hydrogen and helium line ratios. Because 
at wavelengths shorter than 4.05pm the size of the extinction correction becomes very large and 
sensitive, we have not attempted to fit or model the many shorter wavelength line fluxes that 
have been seen. 

The plots in Figure 2 also illustrate that no combination of wind and static atmosphere 
such as we have considered here is able to explain the very strong measured FIR fluxes in the 
H15a, H12a, and H10a lines reported by Strelnitski et al.. Our results confirm that these lines 
do indeed show sign of significant amplification: hence, lasing. 


III. THE VELOCITY COMPONENTS 

The four observed mid-IR lines show velocity structure spanning as much as 170 kms' 1 
(figures la-d). When the lines are fit to a single gaussian profile their peak velocities and their 
widths closely match (see Table 3): V^r * +12±2 kms' 1 , AV(FWHM)=90±8 kms' 1 . This value 
for Vlsr is comfortably close to that of the star: +8kms'' (Thum et al., 1995). The high SNR of 
the data make it tempting to try to attribute some of the velocity structure we observe, especially 
in the H7a line, to the several wind or disk components of the emission. A two component fit to 
the profile of H7a gives +30 kms' 1 (strong component) and -17 kms' 1 (weak component), quite 
close to the values seen in the double-peaked Bry line (Hamann and Simon, 1986), and to the 
weakly defined ~+37kms' 1 feature they find in the Bra line, which they attribute to the surface of 
an inner wind evaporated from the disk. Hamann and Simon show that a bipolar flow viewed 
approximately edge-on, as in the case of MWC 349A, is not by itself able to explain the line 
structure they observe, and conclude some disk-like rotation is needed; indeed this is the origin of 
the submillimeter lines' structures. If we compare the fluxes in the bright (that is, the positive 
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velocity) components of our two mid-IR a-lines with those predicted from our previous wind 
model calculations, we find agreement to within 30%. The high velocity component of the H6a is 
weak, however, and cannot by itself account for the flux from an ionized, static atmosphere 
component (normalized as above to an intrinsic total Bra line flux of 3.0xl0' 17 W/cm 2 ). The 
negative velocity components are themselves weakly subdivided (Table 1): each of the two a- 
lines contains a doublet at high negative velocities separated by about 65±5 kms' 1 . (We note that 
possible contamination from the He I 7-6 line at this velocity should be 10 time weaker.) In a 
Keplerian disk, such as that considered by Strelnitski, Smith and Ponomarev (1966), the velocity 
at 30AU from the star is ~50 kms' 1 (as measured by the H30a maser line velocity and the 
interferometric observations of Planesas et al., 1992). Hence if some of the velocity pair structure 
we see in the IR line emission arises from the disk, this velocity separation suggests a region 
located at about 25-30 AU from the star. Strelnitski, Smith and Ponomarev show that were laser 
emission to be present in the H7a line, it would be expected from a much denser region: at a 
distance of more like 3xl0 12 cm, within the static atmosphere of the Hollenbach model but much 
closer to the stellar photosphere, where the velocity separatin of the red and blue features would 
be much greater than that observed. Alternatively the dense regions might exit in knots of 
material located further out in the disk. Additional observations and temporal monitoring would 
be helpful in sorting out this possibility. 


IV . CONCLUSIONS 

We present new observations of four mid-IR hydrogen recombination lines, whose 
strengths are consistent with conventional emission from the ionized atmosphere around the 
Keplerian disk of the star, the strong bipolar wind, or contributions from both. There is even mild 
agreement in the observed velocity structure of the lines. No amplification (lasing) is required in 
these new lines reported here, and the results confirm and strengthen our earlier conclusions that 
the FIR lines are lasing. The detailed conditions needed in the source to produce detectable laser 
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emission in the mid-IR lines are discussed in Strelnitski, Ponomarev and Smith (1966). Although 
the theory of Hummer and Storey (1995) indicates inversions can be present between the levels of 
all the transitions we report here, considerations of laser saturation and/or possible beaming may 
make it hard to produce or detect significant amplification (Strelnitski, Ponomarev and Smith). 
Alternatively it might just be that the densities needed for inversions at, say the HI a line -- 
between about 2xl0 10 and 2x10" cm' 3 (op cit .) — don't exist as a distinct disk annular region, 
perhaps because at these high values the material is too close to the star to retain its Keplerian 
character (see the "weak wind" model of Hollenbach etal., 1994). These results fill in a portion 
of the critical infrared wavelength interval within which this remarkable source switches from 
conventional emission to shining in laser lines from its disk. 
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Figure Captions: 


Figure 1 a-c: Irshell spectra of the four hydrogen recombination lines detected. Solid and dotted 
lines show the result of gaussian fitting. The locations of weak telluric features are also indicated. 

Figure 2: Observed mid- and far-IR lines, also including the Bra line. The solid line shows the 
predicted flux in the Hna lines from a bipolar wind of constant velocity, normalized to the 
observed flux in the H92a transition at 3.1 cm. The dotted line shows the predicted flux in the 
Ha lines from both the wind plus a component from a thin HH-region, with the latter contribution 
scaled to the flux in Bra. The three FIR lines are from Strelnitski et al., 1966. The x-axis is 
marked with diamonds showing the wavelength locations of all the Hna lines (H16a - H4a); the 
theoretical curves were produced by calculating fluxes for each of these lines, and connecting 
points with straight line segments. 
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Table 1. Log of Observations 


Transition 

A.(nm) 

UT Date (1994) 

'^'“integration («*) 

Resolution (kms 1 ) 

H7a 

19.062 

June 7 

1258 

15 

H6a 

12.372 

June 7 

419 

25 

H6a 

12.372 

June 8 

1258 

25 

H7p 

11.309 

June 8 

210 

15 

3 Total on source integration time 


Transition 

Table 2. Integrated Line Fluxes 
A(|im) Integrated Line Fluxes 3 (x 10 15 Wm' 2 ) 

H7a 

19.062 

8.7 

H6a 

12.372 

22.0 

H7p 

11.309 

6.5 

H8y 

12.387 

3.9 


a Systematic uncertainties are estimated to be ~ 15%; relative uncertainties are <5%. 
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Table 3. Line Component Parameters 


Transition 

^LSR 

(kms' 1 ) 

^^FWHM 

(kms 1 ) 

Flux 

(xlO 15 Wm' 2 ) 

H7a 

30 ±2 

57 ± 3 

4.8 ± .3 


-17 ±2 

48 ± 4 

2.5 ± .3 


-74 ±4 

33 ± 11 

0.2 ±.l 


-144 ±5 

95 ± 13 

0.9 ±.l 

(single gaussian fit: 13.5 

82.4 ± 3 

8.7 ) 

H6a 

13 ± 1 

90 ± 1 

17.4 ± 0.3 


-103 ±7 

35 ± 32 

0.3 ± 0.3 


-160 ± 17 

116 ± 31 

1.9 ± 0.6 

H7p 

14. 3± 3 

82. 6± 4 

6.5 ± 0.6 

H8y 

11 ±2 

98 ±6 

3.9 ± 0.3 
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